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Abstract 


This paper considers the leader-follower tracking control problem for 
linear interconnected systems with undirected topology and linear dy¬ 
namic coupling. Interactions between the systems are treated as linear 
dynamic uncertainty and are described in terms of integral quadratic con¬ 
straints (IQCs). A consensus-type tracking control protocol is proposed 
for each system based on its state relative its neighbors. In addition a se¬ 
lected set of subsystems uses for control their relative states with respect 
to the leader. Two methods are proposed for the design of this control 
protocol. One method uses a coordinate transformation to recast the 
protocol design problem as a decentralized robust control problem for an 
auxiliary interconnected large scale system. Another method is direct, it 
does not employ coordinate transformation; it also allows for more general 
linear uncertain interactions. Using these methods, sufficient conditions 
are obtained which guarantee that the system tracks the leader. These 
conditions guarantee a suboptimal bound on the system consensus and 
tracking performance. The proposed methods are compared using a sim¬ 
ulation example, and their effectiveness is discussed. Also, algorithms are 
proposed for computing suboptimal controllers. 
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1 Introduction 


Theoretical study of distributed coordination and control has received in¬ 
creasing attention i n the past decade, d ue to its broad appl ications in unmanned 
air vehicl es (UAVs) (Beard et al. , 2002 ), formation control ( Fax fc Murravl 2004), 
flocki ng (lOlfati-Saberl . 1200611 and distributed sensor networks (ICortes fc Build . 

l2003h . etc. As a result, muc h progress has been made in the study of cooperative _ 

contr ol of complex systems (lOlfati-Saber. Fax, fc Murravll2007l : lRen. Beard, fc: Atkinsl . 
[20^, with the aim to develop feedback control tools to achieve a desired system 
behavior. In particular, synchronization problems f or interconnected networks 
of complex dynamical systems are actively studied ( Arenas et ah . 20081 : Tuna, 

2008112OOOII . 


There exist a number of approaches to achieve synchronized behavior in 
systems comprised of many dynamic s ubsystems-agents. These approache s in¬ 
clude the average consensus approach ( Olfati-Sa ber. Fax, fc Murray. 2007 |). the 
appro ach based on internal model principle (fWieland. Sepulchre, fc Allgower . 
201 ih . and the leader-follower approach. In the latter approach, one of the 


agents is designated to serve as a leader, and interconnec tions within the system _ 

are d e signed to let the re s t of th e system follow the leader (iGrip. Yang, Saberi. fc Stoorvogel . 


2OI2I : Pecora fc Carroll , 1990ll . This approach, known as leader-follower ap¬ 


proach, is the main focus of this paper. 

The majority of leader-follower problems considered in the literature as¬ 
sume that dynamics of the agents are dynamically decoupl ed ( Hong. Hu. fc Gaj . 
20061 : Jadbabaie. Lin, fc Morsel 120031 : T^ao. Li. fc Duanl. 2013 ), and th e infor¬ 


mation flow between the subsystems is directed (IRen fc Atkinsl . l2007^ and is 
used for control. While these assumptions are justifiable in the case of multi¬ 
agent systems such as autonomous yehicle formations, in many physical sys¬ 
tems, interactions are unayoidable and haye undirected nature (ISiliakl . Il978 ; 


Siliak fc Zeceyid . I2OO5II . The Newtonian interaction between mechanical sys¬ 
tems (e.g., the grayitational attraction between satellites), and the Coulomb 
forces between charged particles are the examples of such undirected inter¬ 
actions. Moreoyer, implementation of control protocols using these physical 
principles (e.g., by interconnecting physical masses with springs and dampers) 
ineyitably leads to undirected control interactions. Hence, there is a need to 
explore situations where undirected interaction s occur at both interconnection 
and control leyel (jPersis. Sailer, fc Wirthl . 120131) . This motiyates us to consider 
undirected control schemes. 

Gompared with the existing work in the field of the leader-follower track¬ 
ing consensus problem, we consider a quite general class of physical inter¬ 
actions between subsystems. These interactions include both static and dy¬ 
namic interactions, such as unmodelled linear dynamics, uncertain input de¬ 
lays and norm-bounded uncertainties. To capture such a broad class of in¬ 
teractions, we regard them as an uncertainty and describe them in terms of 
time-domain inte gral quadratic inequa l ities known as Integral Quadratic Con¬ 


strain ts (IQCs) (jMegretski fc Rantzeii . 119971 : iPetersen. Ugrinoyskii. fc Saykin . 
2OOOII . The IQC modeling is a well established technique to describe uncer- 
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tain interactions between subsystems in a large scale system. It has led to 
a nu mber of solutions to optimal and suboptimal decentralized control prob¬ 

lems (Li. Ugrinovskii. fc OrsiL 2007 : Ugrinovskii et ah . 2000 l : Ugrinovskii fc Pota . 


2005h . 


As in these references, the IQC modelling allows us to account for the effects 
of interconnections between subsystems from a robustness viewpoint. However, 
different from the above references, the IQC methodology is developed here for 

the design of distributed consensus-type feedback tracking con trollers. _ 

_ the context of robust consensus analysis, the recent paper ()Trentelman. Takaba. fc Monshizadeh , 

20I3h is worth mentioning, which considers robust consensus protocols for syn¬ 


chronization of multi-agent systems under additive uncertain perturbations with 
bounded Hoc norm. Since the IQC conditions in our paper capture uncertain 
perturbations with bounded L 2 gain, we note a similarity between the two 
uncertainty classes. However, thanks to the time-domain IQC modelling, our 
paper goes beyond the analysis of robust consensus. It develops the technique 
for leader-follower distributed tracking control synthesis, which provides an op¬ 
timized guarantee of perform ance of the leader-follower tracking system under 
consideration (note also that iTrentelman. Takaba. k, Monshizadeh! ( 201311 con¬ 
sider a leader less network). 

The key element of our approach to the leader-follower tracking control 
synthesis is an optimization formulation which imposes a cost on the worst- 
case consensus tracking performance of the system, as well as on protocol ac¬ 
tions j_This_a££roanhHshis 2 ired_by_rtie_ 2 ecentj;esultson_distributed_JjQR^_de; _ 

sign (Borrelli fc Keviczkvl 2008 : Zhang. Lewis, fc Daj . l20IlllZhao. Duan. Wen fc Cheiil . 


20I2[) . It allows us to recast the original consensus tracking problem as a de¬ 


centralized guaranteed cost control problem for a certain auxiliary large-scale 
system. This leads to a distributed control design method for the system of 
coupled subsystems, where local tuning parameters can be chosen to minimize 
the bound on the consensus performance of the protocol leading to a suboptimal 
guaranteed performance. This reduces the original problem to an optimization 
problem involving coupled parameterized linear matrix inequalities (LMIs). We 
also show that the design of the tracking protocol can be simplified using decou¬ 
pled LMIs. This however leads to a weaker tracking result in that we can only 
guarantee a greater bound on consensus tracking performance. Furthermore, we 
compare this method with an alternative method based on direct overbound¬ 
ing of the original performance cost. The advantage of this method is that it 
can be extended to the case of interconnected systems with more general linear 
uncertain dynamic coupling, as demonstrated in Section [3.41 

The main contributions of the paper are sufficient conditions for the design 
of a guaranteed consensus tracking performance protocol for interconnected sys¬ 
tems subject to linear dynamic IQC-constrained coupling. To derive such condi¬ 
tions, we first transform the underlying guaranteed consensus performance con¬ 
trol problem into a guaranteed cost decentralized robust control for an auxiliary 
large scale system, which is comprised of coupled subsystems. The intercon- 
nections pose an additiona l difSc u lty here, compared w it h rec ent results, e.g.. 


iLi. Duan. Chen, fc Huang ( 2010h : Zhang. Lewis, fc Dad ( 20Ilh . where similar 
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transformations resulted in a set of completely decoupled stabilization problems. 
To overcome the effect of the interconnections, we emp loy the minimax control 
design methodology of decent r alized control synthesis (iLi. Ugrinovskii. fc OrsiL 
200'H Ugrinovskii et ah . 200(1 Ugrinovskii fc Pot^ 2005ll . We then discuss an 


alternative sufficient condition whose derivation does not involve the coordinate 
transformation. We show using an example that our main result may offer an 
advantage, compared with this alternative condition. Finally, the computational 
algorithms are introduced to optimize the proposed guaranteed bounds on the 
consensus tracking performance. 

The preliminary version of the paper was p resented at th e 2013 American 
Contr ol Conference ( Cheng fc Ugrinovskiil l2013f) . Compared to Cheng fc Ugrinovskii 
(201^, this paper has been substantially extended. Firstly, in this paper, more 
general linear uncertain coupling is considered, and the leader is allowed to dy¬ 
namically couple with some of the followers. In addition, we present a detailed 
comparison of the results in ICheng fc Ugrinovskii ( 2013ll with those obtained 
using a direct technique which does not involve coordinate transformation. An¬ 
other extension in this paper is concerned with the computational algorithms, 
which demonstrate how the design of a suboptimal tracking protocol can be 
carried out by minimizing the proposed guaranteed bound on the consensus 
tracking performance. 

The paper is organized as follows. Section[5]includes the problem formulation 
and some preliminaries. The main results are given in Section [31 In Sectional 
the computational algorithms are introduced. Section [5] provides the illustrative 
example. The conclusions are given in Section |6l All the proofs are given in the 
Appendix. 


2 Problem Formulation and Preliminaries 


2.1 Interconnection and communication graphs 

Unlike many pap ers that study the l eader -f ollower tracking prob lem for de¬ 
coupled systems (cf. iHong. Hu. fc Gao ( 20061 1: [Ren fc Atkins! ( 2007ll ). we draw 
a distinction between the network representing ’physical’ interactions (including 
the leader) and the network that realizes ’control’ interactions. The rationale 
for considering the two-network structure is twofold. Firstly, synchronization 
protocols must be designed for followers only, and should have no direct im¬ 
pact on the leader. Also, the control interactions do not have to replicate the 
topology of physical interconnections. 

Consider an undirected interconnection graph Q = where V = 

{0,..., A^} is a finite nonempty node set and £1'^ C V x V is an edge set of 
unordered pairs of nodes. Without loss of generality, the node 0 will be assigned 
to represent the leader, while the nodes from the set Vo = {!,...,A^} will 
represent the followers. The coupling between the followers is described by a 
subgraph Qo of G defined on the node set Vo with the edge set £q C Vo x Vo. 
The edge (z, j) in the edge sets £’^, £q means that nodes i and j influence each 
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other through physical interconnections. 

Let Aff be the adjacency matrix of the subgraph Qo', o,fj are the elements of its 
zth row where afj = 1 if (i,j) G and afj = 0 otherwise. The Laplacian matrix 

of the subgraph Qo is defined as £g = T — A^, where T = diag{/i,..., /a?} G 
j^NxN jg in-degree matrix of Go, i.e., the diagonal matrix, whose diagonal 
elements are the in-degrees of the corresponding nodes of the graph Go, fi = 
for i = 1,N. In accordance with this structure, the adjacency 
matrix A'^ of the undirected graph G is obtained by augmenting Ao as follows 



where d = [di ... d^]', with dj = 1 if there is the interconnection between the 
zth follower and the leader, and = 0 otherwise. 

Also, consider an undirected control graph C = (VoyS^^jA^^) with the same 
vertex set Vo and an undirected edge set £1° C Vo x Vq. An unordered pair (z, j) 
in the edge set £1° indicates that nodes z and j obtain information from each 
other, which they will use for control. C is assumed to have no self-loops or 
repeated edges. The adjacency matrix A^^ = [a^j G jiNxn undirected 

graph C is defined as = 1 if (z, j) G and = oG = 0 otherwise. 

The degree matrix TL = diag{hi,..., Hn} G is a diagonal matrix, whose 

diagonal elements are hi = for z = 1,..., TV. The Laplacian matrix of 

this graph is denoted as = H — A‘^. It is symmetric since C is undirected. 

We assume throughout the paper that the leader is observed by a subset of 
followers. If the leader is observed by follower z, we extend the graph C by adding 
the directed edge (z, 0), and assign this edge with the weighting gt = I, otherwise 
we let (/i = 0. We refer to node z with ^ 0 as a pinned or controlled node. 
Denote the pinning matrix as G = diag [^fi,..., G The system is 

assumed to have at least one follower which can observe the leader, hence G 0. 
The extended graph represents the communication topology for control and is 
denoted as C. Let g =[gi ... gN]', its adjacency matrix A‘^ is dehned as 



Finally, we introduce the notation for neighborhoods in the above graphs. 
Node j is called a neighbor of node z in the graph C {G or Go, respectively) 
if (z, j) G £1° or £q, respectively). The sets of neighbors of node z in the 

graphs C, G and Go are denoted as Sf = {jj(z, j) G £1°}, Sf = {j\{i,j) G E'^}, 
and Si = {j\{i,j) G Eq}, respectively. 

2.2 Problem Formulation 

Consider a system consisting of A^ -I- 1 interconnected subsystems; these 
interconnections are described by the undirected graph G- Dynamics of the zth 
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subsystem are described by the equation 


Xi = Axi + BiUi + ^2 ^ y)(t,Xj(.)|o - a:i(.)lo), (1) 

j^sf 


where the notation ?/(.)|g) describes an operator mapping functions y{s), 
0 < s < t, into Also, Xi G 5R” is the state, Ui G is the control input. 
We note that the last term in © reflects a relative, time-varying nature of 
interactions between the subsystems. 

Let L9.[0 ,oo) be the space of functions y(.) : [0,oo) —>■ 5R" such that 
/g‘||y(t)fdt<oo, Vt>0. 

Assumption 1 Given a matrix C G the mapping <y9(.,.) satisfies the 

following assumptions: 

(i) Vy G L^jO,oo), (p(.,y(.)| 5 ) G L™jO,oo). 

(a) Vt > 0, (p{t,y) is linear in the second argument; i.e., if y = aiyi -I- 022 / 2 , 
then (p{t,y{%) = ai(p{t,yi{.%) + a 2 (p{t,y 2 {.)\o). 

(Hi) There exists a sequence {ti}, U —>■ 00 , such that for every ti, the following 
IQC holds 

[ \\T{t,y{-)\lWdt< [ WCyW^dt, Vy G L2e[0,oo). (2) 

Jo Jo 


The class of such operators will be denoted by 


Remark 1 Assumptionl^ captures some common classes of uncertain coupling. 
For example, (p can be a linear causal operator from the Hardy space Hoc- 
Such ope r ators have extension to operators mapping L2e[0,oo) into L2e[0,oo) 
' Willemli . I 97 A) . For instance, it is easy to show that unmodelled dynamics 


described as 


f Ci = -aiO + y(t), Ci(0)=0, 
I T(t,y(-)io) = Ci(i), 


satisfy (0). Then the term ip(t,Xj{.) — Xi(.)) in (Qp reduces to (p(xj{.) — Xi(.)) 
and can be interpreted as an action based on relative measurements and applied 
through a stationary dynamic channel with memory. Uncertain input delay in 
receiving relative states is also allowed by this assumption, which can be described 
by choosing 


TiUym 


Cy{t - r), t> T, 

0, 0 < t < T, 


where t is uncertain delay parameter. For this ip we have 

rMt,yi.)\l)fdt= r "||Cy(s)fds< r\\Cy{s)fds. 
Jo Jo Jo 
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This implies that uncertain input delay in receiving relative states is allowed by 
Assumption 1. 

Finally, captures norm-bounded uncertain coupling by allowing the un¬ 
certainty of the form j/(.)|q) = A(t)y{t) where A is a time varying matrix 
such that A'(t)A{t) < C'C. 

Since we have designated node 0 to be the leader, the leader is not controlled, 
i.e., uo = 0. On the contrary, all other follower nodes will be controlled to track 
the dynamics of the leader node. In this paper we are concerned with finding a 
control protocol for each follower node i, of the form 

- x^)-{-g^{xo - Xi)}, (3) 

ie-Sf 


where K G is the feedback gain matrix to be found. 

As a me asure of the system perfor mance, we will use the quadratic cost 
function fcf. iBorrelli fc Keviczkv ( 2008ll '). 

{Xj - XiYQjxj - Xi)-\-giixo - XiYQjxo - Xj) + u'^Ru,)dt, 

(4) 


where Q = Q' > 0 and R = R' > Q are given weighting matrices, u denotes 
the vector u = [u'l ... u^]'. The cost function (jd)) penalizes the system inputs. 
It also penalizes the disagreement between subsystems and their neighbors as 
well as the tracking error between the leader and the pinned subsystems which 
observe the leader. 

The problem in this paper is to find a control protocol ([3|) which solves the 
following guaranteed consensus performance tracking problem: 

Problem 1 Under Assumption[ll find a control protocol of the form m such 
that 


sup J"(m) < oo. (5) 

^0 

It will be shown later that ([S]) implies Ci G L2[0, oo) Vi = 1,... ,A, where 
Ci = xq — Xi is the tracking error at node i. Hence, solving Problem [T] will 
guarantee that all followers synchronize to the leader in the L 2 sense. 


2.3 Associated Decentralized Guaranteed Cost Control Prob¬ 
lem 


In this section, we introduce an auxiliary decentralized guaranteed cost 
contr ol problem for an interconnected large scale system. Our app roach fol¬ 
lows ( Li, Duan, Chen, fc Huanj . 20101 : Izhang. Lewis, fc Dasl . l201lll . however 
here it results in a collection of coupled subsystems. 
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From (ED and taking the linearity of the operator ip into account, dynamics 
of the tracking error vectors satisfy the equation 

e, = Aei - Biu,-B2 ^ ei(.)lo) - <p(t,ej(.)lo)) - B 2 efc(.)lo) - ei(.)lo)- 

jeSi fe: dfc = l 

( 6 ) 


Then the closed loop system consisting of the error dynamics (jHD and the 
protocol dSD can be represented as 


e = {lN® A)e + ((£" + G) 0 (BiK))e - ((£^ + B + B) 0 B2)^(t), (7) 


where 0 denotes the Kronecker product, and e = ... e^]', B = diag[di,..., d^], 


B = 

■ di 
di 

d2 ■ 
d2 ■ 

1 

, m = 

ip{t,ei{.%) 

£(£e2(.)|o) 



d2 ■ 

dN 


. £(£eiv(.)|(,) _ 


It was shown in Hong. Hu. fc Gaol ( 2006ll that if the communication graph C 
is connected and at least one agent can observe the leader, then the symmetric 
matrix + G is positive definite, Hence all its eigenvalues are positive. Let 
T G be an orthogonal matrix such that 


T-\C^ + G)T = J = diag [Ai,..., A^]. (8) 


Also, let e = (T ^0 In)^, e = [e'^ ... e^]' and ^'(t) = (T ^ 0 Using 

this coordinate transformation, the system ED can be represented in terms of 
as 


e =(7jv 0A + J0 {BiK))e - {M 0 H 2 )«'(t), (9) 


where M = T ^(£q + B + B)T and 




N 

N 

i=i 


£(££i(-)Io) 


Here we used the assumption that ip{t, •) is a linear operator. It follows from 
([HD that the system ([HD can be regarded as a closed loop system consisting of N 
interconnected linear uncertain subsystems of the following form 


ii — Asi + BiiUi + Ei^i + LiTji, 


( 10 ) 
















each governed by a state feedback controller Ui = Ksi. Here we have used the 
following notation 


( 11 ) 

Vi = [^1 ■ ■ ■ Ci-1 Ci+1 ■ • ■ ^^v]^ (12) 

-Bii = Ai-Bi, 

Li = —B2[Mi^iI, ..., , Mi^ml]- 

From Assumption [U the following two inequalities hold for all i = 1,..., A^: 


rm?dt< r\\Ce,fdt, 

Jo Jo 



(13) 

(14) 


It follows from m and (TH)) that the collection of uncertainty inputs ^i, rji, i = 

1,..., N, represents an admissible local uncertainty and a dmissible interconnec¬ 
tion inputs for the large-scale system (ITOt resp e ctively; see Li. Ugrinovskii. fc Orsil _ 

(20(H); Petersen. Ugrinovskii. &: Savkini ( 200^ : Ugrinovskii et al.l ( 2000ll : Ugrinovskii fc Pota 
(2005j). Let S, H be the sets of all uncertainty inputs and admissible intercon¬ 
nection inputs for the system (fTUl) for which conditions (fT^ and (ITU) hold. 

Thus, we conclude that if (p satisfies the conditions in Assumption [TJ then the 
corresponding signals (fTTI) . (IT^ belong to 5, H, respectively. 

Next, consider the performance cost (jH). It is possible to show that 

J{u)=j^ (e'((£^-f G)®Q)e-f it'(/(g)i?)it)dt. (15) 


Since T is an orthogonal matrix and e = {T ^ (8) /n)e, then e = {T ® In)e 


and 


u=-[{{C‘= + G)T)(E}K]e. 

Since T'T = TT' = In, this allows the performance cost to be expressed as 

J(n) = ^ (e'((T'(/:" + G)T) 0 g)e + + G)iC‘= + G)T) 0 {K'RK))e')dt 

^ POO 

= {X^e',Qei + Xje',K'RKei)dt. (16) 

i=ldO 

Thus we conclude that for u = —((£° -I- G) 0 K)e and u = {u[,... ,u'j^y, 

Ui — KSi, 

J{u) = J{u), (17) 
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where 


^ poo 

J(u) = W/ {Xie'^Qe^ +X'fu^Rui)dt. (18) 

i=i-^o 

Now consider the auxiliary decentralized guaranteed cost control problem 
associated with the uncertain large scale system comprised of the subsystems 
cni), with uncertainty inputs (HU and interconnections subject to the 

IQCs (fT^ . (fHl) . In this problem we wish to find a decentralized state feedback 
controller u = {u [,..., Ui = Ksi such that 


sup < oo. (19) 

H,n 

The connection between this problem and Problem [T] is given in the following 
lemma. 

Lemma 1 Under Assumption Ql if the decentralized state feedbaek controller 
u = {u'l,... ~ solves the auxiliary decentralized guaranteed cost 

control problem for the collection of systems m and the eost funetion m, 
then the control protocol with the gain matrix K solves ProblemUi 

The proof of the Lemma and all other results are given in the Appendix. 

Note that since £“ + G is positive definite, then it follows from (HID and dUD 
that the protocol ([HD with the gain matrix K obtained from the auxiliary de¬ 
centralized control problem will also guarantee e € L 2 [ 0 ,oo). 

Remark 2 The system transformation described in this section reduces the sys¬ 
tem to a collection of interconnected systems m where each node must know 
its corresponding eigenvalue of the matrix CT P G. When the graph topology is 
completely known at each node, these eigenvalues can be readily computed. But 
even if the graph topology is not kno wn at each node, these eigenvalues can be 
estim ated in a decentralized manner t Franceschelli. Gasvarri. Giua. & Seatzv . 

p.ni.i) , 

3 The Main Results 

3.1 Sufficient Conditions for Gnaranteed Performance Leader- 
follower Tracking Control 

The main results of this paper are sufficient conditions under which the 
control protocol ([SD solves the guaranteed consensus performance leader-follower 
tracking control problem. The first such condition is now presented. 
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Theorem 1 If there exist matrices Y = Y' > Q, Y € F € and 

constants tti > 0, 6 i > 0, i = such that the following LMIs (with 

respect to Y, F, ^ and ^) are satisfied simultaneously 



F' 

YQy^ 

YC' 

1' (g) YC ' 


F 

1 n-l 

0 

0 

0 


qy^Y 

0 

-I 

0 

0 

< 0 

CY 

0 

0 

-d-/ 

TTj 

0 


1®CY 

0 

0 

0 

-©z”' . 



( 20 ) 


where 1 = [1 ... 1]' S ^, Qi = >^iQ, ©z = diag[6>i/,..., 0i_i J, 9i+il,6^1], 
and 


.Mt 


T m2. 

2 




Zi — AY + Y A' + XiF' B'l + XiBiF + (—’■—h 


-)B2B'^, 


then the control protocol 0) with K = FY~^ solves Problem\^ Furthermore, 
this protocol guarantees the following bound on the closed loop system perfor¬ 
mance 


N 

sup J{u) < E e'(0)y-ie,(0). (21) 

^0 i=l 

In the special case, when there is no interconnection between the subsystems 
(i.e., i ?2 = 0 and C = 0 in ([T])), with F = —{1/X)R~^B[ and A = maxAi, the 

i 

result of Theorem [1] reduces to the following Corollary. 

Corollary 1 Consider the case B 2 = 0 and C = 0. If there exist matrix 
X — X' > 0, X G such that the following Riccati inequalities are satisfied 

simultaneously 

XA + A'X - ^(2 - ^)XBiR-^B[X + XtQ < 0, (22) 

A A 

then the control protocol (0) with K = —(1/A)i? ^B[X solves Problem Ui for 
the corresponding systems of decoupled subsystems. Furthermore, this protocol 
guarantees the following bound on the closed loop system performance 

N 

supj(w) < ^e'(0)Xe,(0). (23) 

=0 i=i 

Rema rk 3 The proposed condition Y2‘A) is similar to the ARE condition in ref¬ 
erence IZhana. Lewis. & Dak \20li) . The difference between the two conditions 
is due to including the performance specification in our design. 
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3.2 Simplified Sufficient Conditions for Guaranteed Per¬ 
formance Leader-follower Tracking Control 

According to Theorem [I] one has to solve N coupled LMIs to obtain the 
control gain K. To simplify the calculation, it is possible to require only one 
LMI to be feasible, as follows 


Z r(AQ)i/2 YC YC 


(AQ)i/2y _/ 0 0 

CY 0 -i/ 0 


(24) 


0 0 . 


where A = min Ai, = maxMA, = max ^ Aff -, and 



Unlike the LMIs (j20|) . the LMI (1241) is identical for all nodes, it does not involve 
variables from other nodes’ LMIs. This LMI can be solved at each node inde¬ 
pendently. We show in this section that this enables the control protocol to be 
synthesized at each node in a distributed fashion, resulting in the same protocol 
matrix K for all subsystems. First we present the following theorem. 

Theorem 2 Given R = R' > 0 and Q > 0, if the LMI \24^ in variables 
Y = Y' > 0, Tr~^ > 0 and 0~^ > 0 is feasible, then the control protocol m with 
K = -XX-^R-^B[Y-^ solves Problernm Furthermore, the bound m on the 
closed loop system performance holds with Y obtained from 

The tracking protocol ([3]) requires all subsystems to use the same gain K. 
In Theorem [T] a common gain was obtained because the LMIs ((SHI) are coupled. 
In this section, each node has to solve its own version of the LMI (|24)l . which 
are not coupled. Hence, for all nodes to obtain the same gain K, they must 
compute a common matrix Y and constants tt, 6 . This can be done using the 
following consensus algorithm. 

• Let each node i, i = 1,... ,N, solve the LMI to obtain a feasible 

matrix and constants . 

• Then, for a constant /3, 0 < /3 < l/(maxi=i,..._ 7 v ^i)j and k = 0,1,..., 


define 
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( 25 ) 





0 Y, 

3^S2 



Suppose t he control graph C is strong l y con nected. Then according to 


Theorem 2 of lOlfati-Saber. Fax, k, Murravl (l2007l) . if 0 < /3 < 


then 


limfe^oo = F, limfc_j. 
do) 


OO 

(0) 


max(/ii 

= TT, and limfe_j.oo = 0 exist and are equal to 

Y = J2i ^ = Ei and 0 = J2^ where 0 < < 1, X)* = 1- 

Since the feasibility set of the LMI (IMl) is convex, these matrix Y and constants 
TT, 9 are a feasible solution to the LMI (l24l) . We observe that all nodes converge 
to this solution asymptotically, using the consensus algorithm (1251) . Hence, us¬ 
ing this solution, they compute the common gain matrix K with an arbitrary 
accuracy. 


3.3 An Alternative Approach to Derivation of Distributed 
Controller 

The key technique in the previous discussion was the coordinate transforma¬ 
tion, which enabled the synthesis of the leader-follower tracking control for the 
original interconnected system o to be recast as a decentralized robust control 
problem for an auxiliary interconnected large scale system. It is also possible 
to propose an alternative method, which does not involve such a coordinate 
transformation. In this section, we compare the two techniques. 

The derivation of the leader-follower tracking feedback control proposed in 
the previous sections was based on the following upper bound on the cost func¬ 
tion (HI) 

sup J”(u) < sup (u). (26) 

=0 s,n 

There is an alternative ‘direct’ way to obtain a bound on the cost O as follows 


J{u) <X [ e' (In Q + In (XK'RK)) edt = XV] [ 
Jo ^ Jo 


e'(Q-f AiL'i?iL))eidt. 

(27) 


Then we have 


V .OO 

sup 77 ( 14 ) < A sup ^ / e'^{Q + XK'RK))eidt. (28) 

It is important to note that the expression on the right hand-side of (l27)) can 
also be obtained using the coordinate transformation discussed earlier. However, 
in (1281) the supremum of this quantity is taken over a smaller set Sg of operators 
satisfying the IQC condition ([2]). On the contrary, the auxiliary control problem 
used in the proof of Theorem [U involves the supremum over a larger set of 
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uncertainties, described by the IQC conditions (fT^ and (HI ; see (1^ . Thus 
the two techniques can potentially lead to different upper bounds on the leader 
tracking performance. 

In order to formulate the synthesis result based on the alternative upper 
bound on the performance cost, we first introduce the following matrices. When 
follower i is coupled with the leader then di 0. For those followers, consider 
a matrix Y = Y' > 0, and a collection of positive constants fiij, j € Si, viq 
and /ioi, and define a matrix 11^ 



r A, 


YC 

1' (g) YC 

YC 

YC 1 


Q^/'^Y 

-I 

0 

0 

0 

0 

n. = 

CY 

0 

-^I 

0 

0 

0 

li®CY 

0 

0 

-XI, 

0 

0 


CY 

0 

0 

0 


0 


CY 

0 

0 

0 

0 

NiioJ -1 


where li = [1 ... 1]' G Q = {X/X)Q, Hi = diag[^/,j: i G Sj] and 
Z,=AY + YA'-XBiR-^B[ + {^+Y — + — + y —)B 2 S'. 


In the same manner, for the followers that are decoupled from the leader it 
holds that di = 0. In this case, we consider a matrix Y = Y' > 0, a. collection 
of positive constants Vi and /i^, j € Si, and the matrix 11^ defined as 


n, = 


_Zi YQ^/^ YC l'i®YC' 

Q^/'^Y -JO 0 

CY 0 -^J 0 

U®CY 0 0 -fl* 


(30) 


where Zi is modified to be 


Z, = AY+ YA'-XBiR-^B[ + {^+y —+ y —)B2B'^. 


Theorem 3 If there exist a matrix Y = Y' > 0, Y G and constants 

Vi > 0, fXij > 0, i = l,...,iV, and ViQ > 0, ^oi > 0 (for the nodes with 
di 0) such that the following LMIs (with respect to Y, A ^ y ^ y and 
are satisfied simultaneously 


n, <0, i = i,...,iV, (31) 

then the control protocol m with K = —R ^BfY ^ solves Problem [7J Fur¬ 
thermore, this protocol guarantees the following bound on the closed loop system 
performance 

T2 N 

sup J(m) < — y e'(0)r“^ei(0). (32) 

So A . , 
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Remark 4 Compared with the LMIs i2(A) introduced in Theorem]^ the LMIs 
in Theorem n have different dimensions. The LMIs i20\) have fixed dimension 
{2n + p + Nr) x (2n+p + Nr), where p and r are the row dimension of the input 
Ui and the matrix C, respectively. But the dimensions of the LMIs iSl\) depend 
on di and the in-degrees fi of the nodes. For the nodes coupled with the leader, 
the LMI i31\) has the dimension of (2n+(3 + /i)r) x (2n+(3 + /i)r), and for the 
nodes decoupled from the leader, its dimension is (2n+(l+/i)r) x (2n+(l+/i)r). 
Thus these LMIs are generally smaller than the the LMIs i20i) . Therefore from 
a computational viewpoint, Theorem\^may have some numerical advantage over 
Theorem [7J 

We stress again that the upper bound on the worst-case tracking perfor¬ 
mance is obtained in Theorem [3] using the supremum over a smaller uncertainty 
class than in Theorem [1] However the approach in this section uses a more 
conservative bound on the performance cost function; this leads to a conser¬ 
vative gap between the predicted and actual performance. This gap has been 
demonstrated in the example considered in the next section. 

3.4 Further Extensions 

As another distinction between the two approaches discussed in the previous 
subsections, we note that the approach used in Theorem [3] can deal with inter¬ 
connected systems with more general nonidentical uncertain coupling among 
subsystems. Suppose dynamics of the leader and the ith follower are described 
as 


{ io = Axq + B2 Y. V5ofe(t,a;fc(.)|o - a;o(.)lo), 

k: dfc = l 

= Axi BiUi -I B2 Y - Xi{.)\l)B2d^ipM{t,xo{.)\l - Xi{.)\l), 

feSi 

(33) 

where the notations (pij{t,y{.)\Q), ?/(-)lo) describe linear uncertain opera¬ 

tors mapping a function y{s), 0 < s <t into K™. We note that unlike ([T|), these 
operators are not assumed to be identical, therefore the model (1531) allows for 
nonidentical interconnections between subsystem i and its neighbors. 


Assumption 2 Given a matrix Cij € the mapping ipij satisfies condi¬ 

tions (i) and (ii) of Assumption\^and the following IQC condition: There exists 
a sequence {ti},ti —>■ 00 such that for every ti, 


ll<Py (^:2/(-)lo)ll dt < / \\Cijy\\ dt, Vy G L2e[0, 00 ). 


(34) 


Without loss of generality, we assume that the same sequence {t/} can be chosen 
for all operators ipij, e.g., see ( Li. Ugrinovskii. fc OrsiL[200711 . 
operators will be denoted by Si. Obviously, Sq C Si. 


The class of such 
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Problem 2 Find a control protocol of the form @) such that the system FtS\) 
with this protocol satisfies 


sup fl(u) < oo. (35) 

Si 

For node i, introduce the matrices Ci = [Ch^ ... Ci = J', 

where ji, ■ • ■, are the elements of the neighborhood set Si. 

Similar to Theorem[3l we define the following matrices for nodes i with di 0 
and di — 0, respectively. First, when di 0, consider a matrix Y = F' > 0, 
and a collection of positive constants r'y > 0, pij > 0, j € St, viq > 0, p^i > 0, 
then dehne the matrix F^ 



r 

YQ^!'^ 

FC' 

YC' 




Q^I'^Y 

-I 

0 

0 

0 

0 


CiY 

0 

-w^ 

0 

0 

0 

Fi = 

CiY 

0 

0 

-n. 

0 

0 


CmY 

0 

0 

0 

— -d 

0 


Co^Y 

0 

0 

0 

0 

Nfj.0i^ J 


where Wi = diag[^/, j e Si], ^i = diag[^/, j: i G S'j], and 

Z, = AY+ YA'-XBiR-^B[ +( Y {— + —) + —+ V — 

On the contrary, when di = 0, we consider a collection of positive constants 
and pij, j G Si, and the matrix F^ is defined as 




YQ^C 

YC[ 

YC' 

F, = 

qy^y 

-I 

0 

0 


CAY 

0 

-w. 

0 


CiY 

0 

0 

-Q, 


where Zi is revised as 

Z^ = AY+ YA'-~XBiR-^B[ + (Y {— + —)+ Y — 


Theorem 4 If there exist a matrix F = F' > 0, F G and constants 

Vij > 0, pij > 0, i = and I'm > 0, pen > 0 (when di 0) such 

that the following LMIs (with respect to Y, and Y) satisfied 

simultaneously 

F, <0, i = l,...,N, (36) 
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then the control protocol m with K = —R ^B'lY ^ solves Problem m Fur¬ 
thermore, this protocol guarantees the following bound on the closed loop system 
performance 


\2 ^ 

sup J(m) < — ^e'(0)r"^ei(0). (37) 

— i=l 

The proof of Theorem 0] is similar to the proof of Theorem [3] and is omitted for 
brevity. 


4 The Computational Algorithm 


In this section, we discuss numerical calculation of a suboptimal control 
gain K. According to Theorem [TJ the upper bound on consensus tracking 
performance is given by the right hand side of (EJ. Hence, one can achieve a 
suboptimal guaranteed performance by optimizing this upper bound over the 
feasibility set of the LMIs (1301) : 




{F,F,7ri,e, 


inf 


N 

E 


(l20l) holds} 


e'(0)F-ie,(0). 


(38) 


As in lLi. Ugrinovskii. fc Orsil (l2007l) . the optimization problem 
to minimizing 7 subject to the LMI constraint 


is equivalent 


7 e'( 0 ) 

e( 0 ) In 


> 0 , 


(39) 


where e(0) = [ei(0)' 62 ( 0 )' ... eAr(O)']'- This leads us to introduce the following 
optimization problem in the variables 7 , Y, F, and 6 i. Find 

'^feo]),® 

where the infimum is with respect to 7 , Y, F, iTi and 9i, i = 1,..., N, subject 
to (I 20 )) and (l39)) . We now show that the optimization problem (l38)) and the 
optimization problem (HOI) are equivalent. 


Theorem 5 J*^ - 


In a similar fashion, one can show that the value of the optimization problem 

N 


inf y e'( 0 )Y-ie,( 0 ) 

(gU holds} ^ 


(41) 


is equal to = inf 7 , where the infimum is taken over the feasibility 

set of the LMI (IMl) and (1301) . 
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Theorem 6 J 




Note that it follows from Theorems [H [S] and H] that 

Also, one can show that the value of the optimization problem 


T9 N 


A 


inf :_\'e'( 0 )y 

(ED) holds} A 


-1 


ez(0) (42) 


is equal to the value of the problem Jp 


(ED,® 


= inf 7 subject to ED and 


7 e'(q) 

e( 0 ) In ® (AA-2y) 


> 0 . 


(43) 


Theorem 7 

Based on this discussion, we propose three algorithms for the design of sub- 
optimal protocols of the form EP- The first algorithm is based on Theorems [D 
andEJ 

• Solve the optimization problem (I40L to a desired accuracy, obtaining a 
collection Y, F, iTi, Oi and 7 . It follows from the proof of Theorem [S] that 
(y,F, 7ri,0i) belongs to the feasibility set of the LMIs (|20)) . 

• Using the found Y,F, construct the gain matrix K to be used in (EP, 
by letting K = FY~^. Also, the guaranteed bound on the consensus 
performance of this protocol can be computed, using the expression on 
the right-hand side of equation m- 

The second algorithm follows the same steps, with the exception that the 
first step employs the optimization problem = inf 7 and LMI ESP, 

and the second step uses the value for K given in Theorem E] We present this 
algorithm as a benchmark for Theorems [D and El The third algorithm also 
follows similar steps but uses the optimization problem ([ 4 ^ ~ 7 

LMI 63). 

In each optimization problem considered above the initial conditions of the 
leader and followers are assumed to be known. In practice, the initial states 
of the subsystems may not be known. To circumvent this issue, random initial 
conditions can be assumed to tune the algorithms as was done, for example, in 
( Li. IJgrinovskii. fc Orsill^OOTh . Suppose the initial states of the error dynamics 
are random and satisfy E[ei(0)ei(0)'] = A4, where E is the expectation operator, 
then we have 


E[^e'(0)y-ie,(0) 

2=1 



(44) 
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Tr(.) is the trace of a matrix. Then, taking the first algorithm for example, 
instead of solving the optimization problem (I38L the following optimization 
problem in the variables Y, F,TTi,6i,i = 1... ,N 

minTr(^y-^7w) (45) 

subject to (|20)) can be solved to obtain a control protocol. The second and third 
algorithms can be modified in a similar fashion, when the initial conditions are 
not available. 


5 Example 

To illustrate the proposed design methods, consider a system consisting of 
21 identical pendulums coupled by identical spring-damper systems. Each pen¬ 
dulum is subject to an input as shown in Fig. [TJ Without loss of generality, the 
pendulum labeled 0 is chosen to be the leader and the remaining pendulums are 
the followers. The dynamics of the coupled system are governed by the following 
equations 


mPao = 

mPcii = 

mPa20 = 


— kia‘^{t){ao — ai) — fcia^(t)(ao — 020) — fc 2 a^(t)(do — di) 

— fc2a^(t)(do — d2o) — mglao, 

-kia^{t){ai - Oi-i) - kia'^{t){ai - Oi+i) - k 2 a^{t){ai - di_i) 

—k2a?(t){ai — di+i) — mglai — Ui, i = 1,..., 19, 

-fcia^(t)(Q;20 - aig) - kia'^{t){a2o - ao) - fea^(t)(d20 - dig) 
-fc 2 a^(t)(d 2 o - do) - mgla 2 Q - U 20 , 

(46) 


where I is the length of the pendulums, a{t) is the position of the spring-damper 
along the pendulums, g is the gravitational acceleration constant, m is the mass 
of each pendulum, ki is the spring constant, and fcg is the damping coefficient. 
The position of the spring-damper system can change along the full length of 
the pendulums and is considered to be uncertain, that is 0 < a(t) < Z. 

Choosing the state vectors Xi = i = 0,..., 20, the equation (l46ll can 

be written in the form of o, where 


0 1 
-f oj’ 


Bi 



B2 = 


0 

m 


and — Xi) = ^^-j^[ki k 2 ]{xj — xt). 

Let 5{t) = ^^ 75 ^, C = [ki k 2 ], then (p(t, Xj—Xi) = d{t)C{xj—Xi), 0 < S{t) < 1, 
and the operator ip{t,y) = S{t)Cy satisfies Assumption 1. 

The communication topology of the interconnected system is shown in Fig.[5J 
Note that the subgraph excluding the leader node 0 is undirected. According 
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Figure 1: Interconnected pendulums. 


to this graph, the leader’s position and velocity are available to pendulums 1, 
7, 12 and 18, but are not available to other nodes. Also, all subsystems in this 
example are coupled according to the undirected graph shown in Fig. [3l 

Three simulations were implemented to illustrate the protocol designs based 
on Theorems[l][5]and|31 respectively. We used the same initial conditions for the 
corresponding pendulums in all three simulations and used the same matrices 
Q = [1 0;0 0.1] and R = 0.01. The parameters of the coupled pendulum 
system were chosen to be to = 1kg, I = Im, g = 9.8m/s'^, ki = 0.5N/m, 
k 2 = 0.5fV/ {m/s) and a{t) = 0.5 + 0.4sin(t). 

First, consider the computational algorithm based on Theorems [T] and [SJ 
The problem (l40l) was found to be feasible and yielded the gain matrix K — 
[23.85,40.05]. The simulated relative positions and relative velocities, with re¬ 
spect to the leader, of all pendulums controlled by this control protocol are 
shown in Fig. |4l 

The second simulation and third simulation are based on Theorems [5] and [HI 
and Theorems [3] and [71 respectively by using the same matrices Q and R, and 
the same initial conditions. The control gain matrix K was computed to be 
[205.12,303.52] and [22.72, 77.41], respectively. The simulation results are shown 
in Fig. [Hand Fig. |6l 

Also, for each controller obtained by means of the proposed computational 
algorithms, we directly computed the cost function (jj]). These values are com¬ 
pared with the theoretically predicted bounds on the tracking performance and 
are shown in Table [TJ From the simulation results obtained, compared with 
the method based on Theorem IH the method based on Theorem [3] has much 
larger values of both the theoretically predicted bound on performance and 
the computed performance. This shows that the method based on Theorem [T] 
has a superior guaranteed consensus performance despite a potentially larger 
uncertainty class used in the derivation of the upper bound on the tracking 
performance. Also, the method based on the simplified LMIs of Theorem [2] 
has substantially larger theoretically predicted bound on tracking performance 
compared with Theorems [T] and [31 The computed performance is also inferior in 
this case. Compared with the method based on Theorems [3] and [31 the method 
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Figure 2: Communication graph. 

<!>— 

®~(k) 

Figure 3: Undirected coupling graph. 



based on Theorem [T] enables the followers to synchronize to the leader in a much 
shorter time, with a better guaranteed performance. It is also interesting to note 
that a superior performance in Theorem 1 was achieved using much smaller gain 
values. 


Table 1: Predicted and computed performance of the proposed controllers for 
the uncertain parameter a. 



Control 

Gain 

Predicted 

Bound 

Computed 

Performance 

Theorem [T] 

[ 23.85 40.05] 

19.68 

8.74 

Theorem [5] 

[205.12 303.52] 

3924.87 

341.97 

Theorem [3] 

[22.72 77.41] 

2401.13 

16.46 


6 Conclusions 

Two approaches to the leader-follower tracking control problem with guar¬ 
anteed consensus tracking performance have been discussed in this paper. First, 
the problem was transformed into a decentralized control problem for a system, 
in which the interactions between subsystems satisfy integral quadratic con¬ 
straints. This has allowed us to develop a procedure and sufficient conditions 
for the synthesis of a tracking consensus protocol for the original system. As this 
approach results in coupled LMIs which need to be solved simultaneously, we 
have also presented a result which does not involve coupled LMIs. Furthermore, 
an alternative method has been proposed which does not employ such a trans¬ 
formation, and instead uses overbounding of the performance cost. The latter 
method is shown to allow for an extension to encompass more general intercon- 
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Figure 4: Relative angles (the top figure) and relative velocities of the pendulums 
with respect to the leader, obtained using the algorithm based on Theorems [I] 
and[5l 




Figure 5: Relative angles (the top figure) and relative velocities of the pendulums 
with respect to the leader, obtained using the algorithm based on Theorems!^ 
andlH 
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Figure 6: Relative angles (the top figure) and relative velocities of the pendulums 
with respect to the leader, obtained using the algorithm based on Theorems [3] 
and[3 


nected systems with nonidentical linear uncertain coupling operators. Also, this 
method can be extended to consider th e interconnected systems w i th dir ected 
communication and interaction graphs ( Cheng. Ugrinovskii fc Wenl . 2013 1. 

These design techniques have been compared using an example. It is worth 
noting that the gaps between the predicted performance and computed per¬ 
formance among the three simulation results are considerably different. The 
method based on Theorem [T] exhibits the smallest gap out of the three results. 
The conservatism of Theorem [3] owes to the conservative upper bound on the 
original cost function in (12811 and the particular form of the controller which 
made the inequality in (15711 possible. These upper bounds have shown a no¬ 
ticeable effect in the example. The method based on Theorem [2] appears to be 
significantly more conservative than the methods based on Theorems [T] and [31 
However, Theorem [2] enables the controller gain to be computed in a distributed 
manner, at the expense of degraded performance. 
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8 Appendix 

8.1 Proof of Lemma [1] 

Since the decentralized state feedback controller Ui = Ksi solves the aux¬ 
iliary decentralized guaranteed cost control problem for the collection of the 
systems m, then there exist a constant c > 0 such that 

sup>7({t)<c. (47) 

s,n 

Also we noted that every signal ip which satisfies Assumption [T] gives rise to an 
admissible uncertainty for the large-scale system consisting of subsystems (uni). 
This implies that for any (p G sq, with u = —((/!“ + G) (Si K)e, we have 

J{u) = J{u) < sup J'(u) < c. (48) 

H,n 

Therefore, one obtains sup >/(«) < supj7('u) < c. It implies that the control 
Ho H,n 

protocol ([3]) with the same gain K solves Problem [T] 


8.2 Proof of Theorem [T] 

Using the Schur complement and substituting F = KY, the LMIs (l20l) can 
be transformed into the following Riccati inequality 


Y-^{A + KBiK) -H (A \BiK)'Y-^ + XiK'RK 


Y 


r Mf i jpi 




— + 


L TTi 


Y-^B2B'^Y-^ + [Q, + {it, + 9,)C'C] < 0, (49) 


where 9i = 'Y^9i. 

Consider the following Lyapunov function candidate for the interconnected 
system m- 


N 


V{e) = Y,e[Y- 


£i. 


(50) 


i=l 


For the controller Ui = Ksi, using the Riccati inequality (H31) . we have 


dVje) 

dt 


< - 


N N j^2 

^ - 2e\Y-^M,,,B2^, 




2=1 




2 JA* 


£{Y-^B2B2Y-^£. 


-1, 


6r 


+ 2£'r-iL,ry, - 9, || r?, f +0, || ry, f -9, II Cs, f ). 


(51) 
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By completing the squares on the right hand side of m and using the 
identity 


N N 

i=l j^i 


(52) 


Z = 1 


we obtain 


fti ^ ^ pti 1 , 

/ —^dt < - y / e'^iXfK'RK + Q,)e,dt - V / || f dt 

Jo dt TY! Jo Jo 


^ rti 


i=l 

N 


+ E TT, / ' (II 6 f - II Ce. \\^)dt - E E / " II II' 

i=i do i=lj^^do 

+ E^* / (IIf “EII11^)'^^- 

i=l "'0 j/i 

Here ti is an element of the sequence {<;} from Assumption [T] Finally, using 
the IQCs (fT51) and m and noting that V{e{t)) > 0, we obtain 

N rti 

E / e',{X}K'RK + Q,)edt < H(e(0)). 
i=id^ 

The expression on the right hand side of the above inequality is independent 
of ti. Letting t; —^ oo leads to j{u) < F(e(0)). This conclusion holds for an 
arbitrary collection of inputs rji that satisfy (HU), (HI, respectively. Then 

supj7'(u) < e'(0)y“^ei(0). The claim of the theorem now follows from Lemma[T] 
s,n 
and 


dt 


8.3 Proof of Theorem [2] 

Using Schur complement, the LMI (l24l) is equivalent to the following Riccati 
inequality 


\2 

AY + YA' - ^BiR-^B[ 


w'^ q^' 

LV + yJ 


HzH' + Y[XQ + (tt + 9)C'C]Y < 0, 

(53) 


where 0 = {N — 1)0. 

Since A < Ai < A, substituting F = —■^R~^B[, m = ir, 0i = 0, 0i = (N —1)0 = 
di, Mf. < w"^ and ^ Mf - < q^, then we obtain 


j¥=i 




Y-\A + X,BiK) + (A + X,BiK)'Y-^ + XjK'RK 


+ 


r MA 




— + 


L TTi 


r-iRsH'y-i + [Q, + {tt, + 0^)C'C] < 0. (54) 
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We obtain the Riccati inequality which is equivalent to (1^ . The proof 
then readily follows from Theorem [TJ 


8.4 Proof of Theorem [3] 

The proof is similar to the proof of Theorem [T1 therefore we only present 
the details which are different from that proof. 

When di ^ 0 and Hi is as defined in ([29]), using Schur complement and 
substituting K = —RB[Y~^, the LMIs (IHTl) are equivalent to the following 
Riccati inequality 

Y-\A + \BiK) + {A +XBiKYY-^ + XK'RK 

r.2 _ 1 1 _ 1 

+ (-K y -1-h y - )Y ^B 2 B 2 Y ^ + <3 + (r'i + /ii + Uio + NnQi)C'C < 0, 


where fXi = d-Y- 

3 ■■ ieSj 

When di = 0 and 11^ is defined in dSOD, a similar transformation results in 
the inequality 

Y-\A + XBiK) + (4 + XBiKYY-'^ + XK'RK 

+ (—+y—+ y —)Y-^B2B'2Y-^YQ + {u,Yjii)C'C<Q. (56) 

^ T d-Ok 

jGSi k: dk — 1 

Consider the quadratic Lyapunov function candidate V{e) = 

for the interconnected system comprised of the subsystems ©■ Since K = 

—R~^B[Y~^, we have 


N 


(e. 

i=i jesf 


ej) + gxei) = -2e'((£" + G) ® {Y-^BiR-^B[Y-^Y)e 

N 

<-2'^e',XY-^BiR-^B[Y-^e,. (57) 

i=l 


It follows from (1571) that 


dV{e) 

dt 


N 

< y 2e'y“i (4 + XBiK^ ei 
2=1 


N N 

2y/,e'y-ii? 27 >(i,e.(.)|^) +2y ^ e'r-ii? 2 ^(t, e,(.)|(,) 

2=1 2 = 1 j^Si 


N N 

- 2 E E e'F-ii?2V5(Cefc(.)|(,)-2yd.e'y-iR27>(i,e.(.)|*). (58) 

2=1 k : dk — 1 2=1 
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Then, using the Riccati inequalities (1551) and (1551) . inequality (1551) . and the 
identities 


N N 

i—1 jGSi i—1 

N ^ ^o*e-C"C7ei = N ^ ^ofeCfeC'Cefc, 


i: di — 1 


k: dk — 1 


N 


2—1 i: di — 1 

in a manner similar to the proof of Theorem [1] we obtain the following bound 
Y e',[}.K'RK+ Qy,dt<V{e{G)). (59) 

i=l -^0 


Condition (|32l) now follows from (|28]) since Q = {\/\)Q. 

It also implies that the control protocol (|31) with K = —R~^B'^Y~^ solves 
Problem [TJ 


8.5 Proof of Theorem [5] 

Suppose the LMIs (1211]) and (|39ll have a feasible solution Y,F,TTi,9i and 7 , 
i = 1,..., iV. Then it follows from (l39ll that 

N 

^e'( 0 )y-ie,( 0 ) < 7 . ( 60 ) 

i=l 

Since the feasibility set of the LMIs (1^ . (155)) is a subset of the feasibility set 
of the LMIs ((20l) . then it follows from (l60l) that (135]) ' 

Conversely, for any sufficiently small p > 0, there exist Y,F,0i and tt^, i = 
I,..., verifying (I5TO such that 

N 

( 61 ) 

Let 7 = (j + ei(0)y“^ei(0), where cr > 0 is an arbitrary constant. Then 

7 , Y, F, 9i and tti satisfy conditions (EOj) and (139|). Furthermore < 

7 = cr + Sfci e((0)y“^ei(0) < cr + p + Letting cr, p -)> 0, we have 

(I 35 |) - ■ "^^is completes the proof. 
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